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ABSTRACT 

By combining high-resolution (HST-WFPC2) and wide-field ground based 
(2.2m ESO-WFI) and space (GALEX) observations, we have collected a multi- 
wavelength photometric data base (ranging from the far UV to the near infrared) 
of the galactic globular cluster NGC1904 (M79). The sample covers the entire 
cluster extension, from the very central regions up to the tidal radius. In the 
present paper such a data set is used to study the BSS population and its radial 
distribution. A total number of 39 bright (m.218 < 19.5) BSS has been detected, 
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and they have been found to be highly segregated in the cluster core. No signifi- 
cant upturn in the BSS frequency has been observed in the outskirts of NGC 1904, 
in contrast to other clusters (M 3, 47 Tuc, NGC 6752, M 5) studied with the 
same technique. Such evidences, coupled with the large radius of avoidance es- 
timated for NGC 1904 (r avoi d ~ 30 core radii), indicate that the vast majority of 
the cluster heavy stars (binaries) has already sunk to the core. Accordingly, ex- 
tensive dynamical simulations suggest that BSS formed by mass transfer activity 
in primordial binaries evolving in isolation in the cluster outskirts represent only 
a negligible (0-10%) fraction of the overall population. 

Subject headings: Globular clusters: individual (NGC1904); stars: evolution - 
binaries: close - blue stragglers 



INTRODUCTION 



Blue straggler stars (BSS) appear brighter and bluer than the Turn-Off (TO) point along 
an extension of the Main Sequence in color-magnitude diagrams (CMDs) of stellar popula- 
tions. Hence, they mimic a young stellar population, with masses l arger than the normal 



clust er stars (this is also confirmed by direct mass measurements; e.g. IShara. Saffer fc Livio 



19971 ). BSS are thought to be objects that have increased their initial mass duri ng their 



evolu tion, and two main scenarios have been proposed for their formation (e.g., iBailyn 



19951 ): the collisional scenario suggests that BSS are the end-products of stellar mergers 



induced by collisions (COL-BSS), while in the mass-transfer scenario BSS form by the 



to the complete coalescence of the two stars ( 


Vlateo et al. 


1990 


Pritchet & GlasDev 1 


1991; 


Bailvn & Pinsonneault Ill995: 


Carnev 


Latham 


Tian et al. 


2006; 


Leieh. Sills & Knieee 


2007) 



Hence, understanding the origin of BSS in stellar clusters provides valuable insight both on 
the binary evolution processes and on the effects of dynamical interactions on the (otherwise 
normal) stellar evolution. The MT formation scenario has by recently received further sup- 
port by high-resolution spectroscopic observations, which detected a nomalous Carbon an d 
Oxygen abundances on the surface of a number of BSS in 47 Tuc (IFerraro et al.l l2006al ). 
However the role and relative importance of the two mechanisms are still largely unknown. 



1 Based on observations with the NASA/ESA HST, obtained at the Space Telescope Science Institute, 
which is operated by AURA, Inc., under NASA contract NAS5-26555. Also based on GALEX observations 
(program GI-056) and WFI observations collected at the European Southern Observatory, La Silla, Chile, 
within the observing programs 62.L-0354 and 64.L-0439. 



- 3- 



To clarify the BSS formation and evolution processes we studying the BSS radial distri- 
bution over the entire cluster extension in a number of galac tic globula r clusters (GCs). W e 



completed such studies in 5 GCs: M 3 (IFerraro et al.lll9 97l. 47 Tuc (IFerraro et al.l 120 04). 



NGC 6752 JSabbi et al.ll2004h . u Cen ( IFerraro et al.ll2006bh . and M 5 (ILanzoni et al. 



2007 



see also Warren, Sandquist & Bolte 2006). Apart from a; Cen where mass segregation pro- 
cesses have not yet played a major role in altering the initial BSS distribution, the BSS 
are always highly concentrated in the cluster central regions. Moreover, in M 3, 47 Tuc, 
NGC 6752, and M 5 the BSS fraction decreases at intermediate radii and rises again in the 
outskirsts of the clusters, yielding a bimodal distribution. Preliminary evi d ences of such a 
bimodality have been f ound also in M 55 by IZaggia. Piotto fc Capacciofil (119971 ). Recent 



dynamical simulations (IMapelli et al.l I2004J . 120061 ; ILanzoni et al. 1 120071 ) have been used to 
interpret the observed trends and have shown that a significant fraction ( ^ 50%) of COL- 
BSS is required to account for the observed BSS central peaks. In addition, a fraction of 
20-40% MT-BSS is needed to reproduce the outer increase observed in these clusters. The 
case of uj Cen is reproduced by assuming that the BSS population in this cluster is com- 
posed entirely of MT-BSS. These results demonstrate that detailed studies of the BSS radial 
distribution within GCs are very powerful tools for better understanding the BSS formation 
channels and for probing the complex interplay between dynamics and stellar evolution in 
dense stellar systems. 

In this paper we present multi-wavelength observations of NGC 1904. These observa- 
tions are part of a coordinated project aimed at properly characterize the UV excess of old 
stellar aggregates as globular clusters, in terms of their hot stellar populations, like Hori- 
zontal Branch (HB) and Extreme HB stars, post-Asymptotic Giant Branch stars, BSS, etc. 



From integrated light measurements obtained with UIT (see iDorman. O'Connell fc Rood 



19951 ). NGC 1904 was known to be relatively bright in the UV, and it was selected as a prime 
target in both our high-resolution (using HST) and wide-field (using GALEX) UV surveys. 
We have obtained a large set of data: (i) high-resolution ultraviolet (UV) and optical images 
of the cluster center have been secured with the WFPC2 on board HST; (ii) complementary 
wide-field observations covering the entire cluster extension have been obtained in the UV 
and optical bands by using the far- and near-UV detectors on board the Galaxy Evolution 
Explorer (GALEX) satellite and with ESO-WFI mounted at the 2.2 m ESO telescope, re- 
spectively. The combination of these datasets allowed a study of the structural properties 
of NGC 1904 (thus leading to an accurate redetermination of the center of gravity and the 
surface density profile), and of the radial distribution of the evolved stellar populations (in 
particular the BSS and horizontal branch star distributions have been derived over the en- 
tire cluster extension). While a companion paper (Schiavon et al. 2007, in preparation) will 
focus on the morphology and the structure of the HB, the present paper is devoted to the 
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BSS population. 



OBSERVATIONS AND DATA ANALYSIS 



2.1. The data sets 



The present study is based on a combination of different photometric data sets: 

1. The high-resolution set - It consists of a series of UV, near UV and optical images 
of the cluster center obtained with HST-WFPC2 with two different pointings. In both cases 
the Planetary Camera (PC, the highest resolution instrument with C/046 pixel -1 ) has been 
pointed approximately on the cluster center to efficiently resolve the stars in the highly 
crowded central regions; the three Wide Field Cameras (WFC with resolution O'.'l pixel -1 ) 
have been used to sample the surrounding regions. Observations in Pointing A (Prop. 6607, 
P.I. Ferraro) have been performed through filters F160BW (far-UV), F336W (approximately 
an U filter) and F555W (V), for a total exposure time t exp = 3300, 4400, and 300 sec, respec- 
tively. Pointing B is a set of public HST-WFPC2 observations (Prop. 6095, P.I. Djorgovski) 
obtained through filters F218W (mid-UV), F439W (B) and F555W (V). Because of the dif- 
ferent orientations of the four cameras, this data set is complementary to the former (with 
the PC field of view in common), thus offering full coverage of the innermost regions of 
the cluster (see Figured]). The combined photometric sample is ideal for efficiently studying 
both the hot stellar populations (as the BSS and the HB stars) and the cool red giant branch 
(RGB) population, and to guarantee a proper combination with the wide-field data set (see 
below). 

The photometric reduction o f both the high-resolution sets was carried out using RO- 
MAFOT (IBuonanno et al. Ill983l ). a package developed to perform accurate photometry in 
crowde d fields and specifically optim ized to handle under-sampled Point Spread Functions 
(PSFs; IBuonanno fc Iannicola 1119891). as in t he case of the HST-WFC chips. The standard 
procedure described in iFerraro et al.l (119971 . 2001) was adopted to derive the instrumen- 
tal magnitudes and to calibrate them to the STMAG system by using the zero-points of 



Holtzman et al. I (119951 ). The magnitude lists were finally cross-correlated in order to obtain 



a combined catalog. 

2. The wide-field set - A complementary set of wide-field U, B, and I images was 
secured by using the Wide Field Imager (WFI) at the 2.2m ESO-MPI telescope, during 
an observing run in January 1999 (Progr. ID: 062.L-0354, PI: Ferraro). A set of WFI V 
images (Progr. ID: 064.L-0255) was also retrieved from the ESO-STECF Science Archive. 
Additional deep wide-field images were obtained in the UV band with the satellite GALEX 
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(GI-056, P.I. Schiavon) through the FUV (1350-1750 A) and NUV (1750-2800 A) detectors. 
With a global field of view (FoV) of 34' x 34', the WFI observations cover the entire cluster 
extension. There is also full coverage of the cluster in the UV thanks to the large GALEX 
FoV, which is approximately 1 deg in diameter and includes the WFI FoV (see Figured where 
the cluster is roughly centered on WFI CCD #2). However, because of the low resolution of 
the instrument (4" and 6" in the FUV and NUV channels, respectively), GALEX data have 
been used to sample only the external cluster regions not covered by HST. 

The raw WFI images were corrected for bias and flat field, and the overscan regions were 
trimmed using IRAF@ tools (mscred package). Standard crowded field photometry, including 
PSF modeling, was carried out independently on each image using DAOPHOTII/ALLSTAR 
( IStetson I Il987l ). For each WFI chip a catalog listing the instrumental U, B, V, and I 
magnitudes was obta ined by cross-c o rrelat i ng the single-b and c atalogs. Seve r al hu ndred 
stars in common with iKravtsov et al.l (119971 ). IStetson I (120001 ). and lFerraro et al.l (119921 ) have 
been used to transform the instrumental U, B, V, and / magnitudes to the Johnson/ Cousins 
photometric system. 

As for the WFI data, also for GALEX observations standard photometry and PSF fitting 
were performed independently on each image using DAOPHOTII/ALLSTAR. A combined 
FUV-NUV catalog was then obtained by cross- correlating the single-band catalogs. 



2.2. Astrometry and homogenization of the catalogs 



Ferraro et al. 


( 


2001. 


2003) 



stars in the WFI FoV, and a cross-correlation tool specifically developed at the Bologna 
Observatory (Montegriffo et al. 2003, private communication) has been employed to obtain 
an astrometric solution for each WFI chip. Several hundred GSC-II reference stars were 
found in each chip, thus allowing an accurate absolute positioning of the stars. Then, we 
used more than 3000 and 1500 bright WFI stars in common with the HST and GALEX 
samples, respectively, as secondary astrometric standards, so as to place all the catalogs 
on the same absolute astrometric system. We estimate that the global uncertainties in the 



2 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 



3 Available at 



http : //www-gsss . stsci . edu/Catalogs/GSC/GSC2/GSC2 . htm 
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astrometric solution is of the order of ~ 0"2, both in right ascension (a) and declination (8). 

Once placed on the same coordinate system, the catalogs have been cross-correlated and 
the stars in common have been used to transform all the magnitudes in the same photometric 
system. In particular, the HST STMAG magnitudes have been converted to the WFI ones 
by using the stars in common between the two samples in the optical bands. Then, the 
GALEX FUV and NUV instrumental magnitudes have been calibrated onto the HST itiiqq 
and m2i8 magnitudes, respectively using the stars in common between the GALEX and HST 
samples. 

At the end of the procedure a homogeneous master catalog of magnitudes and absolute 
coordinates of all the stars included in the HST, WFI, and GALEX samples was finally 
produced. 



2.3. Center of gravity and definition of the samples 

Once the absolute positions of individual stars have been obtained, the center of gravity 
Cgrav of NGC 1904 has been determined by averaging the coordi nates a and 8 of all s tars ly ing 



in the PC FoV, following the iterative procedure described in iMontegriffo et al. I (119951 . see 
also Ferraro et al. 2003, 2004). In order to correct for spurious effects due to incompleteness 
in the very inner regions of the cluster, we considered two samples with different limiting 
magnitudes (V < 19 and V < 20), and we computed the barycenter of stars for each 
sample. The two estimates agree within ~ 1", setting C gvwj at a(J2000) = 05 h 24 m ir?09, 
<5(J2000) = — 24° 31' 29"00. The newly determined center of g ravity is locate d at ~ 7" south- 



est (Aa = 7/3, AS = —2") from that previously derived by lHarris I (119961 ) on the basis of 
the surface brightness distribution. 

In order to reduce spurious effects in the most crowded regions of the cluster due to the 
low resolution of the WFI and GALEX observations, we considered only the HST data for 
the inner 85" from the center, this value being imposed by the geometry of the combined 
WFPC2 FoVs (see Figured]). Thus, in the following we define as HST sample the ensemble 
of all the stars observed with HST at r < 85" from C grav , and as External sample all the 
stars detected with WFI and/or GALEX at r > 85", out to ~ 1100" (see Figure EJ. The 
CMDs of the HST and External samples in the (V, B — V) planes are shown in Figure [3J 

Note that only the data suitable for the study of the BSS population will be considered 
in the following, while those obtained through filters F160BW and FUV on board HST and 
GALEX, respectively, will be used in a forthcoming paper specifically devoted to the analysis 
the HB properties (Schiavon et al. 2007). 
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2.4. Density profile 



Considering all the stars brighter than V = 20 in the combined HST+External catalog 
(see Figure [3]), we have determined the projected density profile of NGC 1904 by direct 
star counts oyer the entire cluster extension. Following the procedure already described in 
Ferraro et al.l (Il999al . 2004), we have divided the entire sample in 31 concentric annuli, each 



centered on C grav and split in an adequate number of sub-sectors (quadrants for the annuli 
totally sampled by the observations, octants elsewhere). The number of stars lying within 
each sub-sector was counted, and the star density was obtained by dividing these values by 
the corresponding sub-sector areas. The stellar density in each annulus was then obtained 
as the average of the sub-sector densities, and the standard deviation was estimated from 
the variance among the sub-sectors. 

The radial density profile thus derived is plotted in Figure HI and the average of the 
three outermost (r > 8'3) surface density measures has been adopted as the background 
contribution (corresponding to 0.95 arcmin -2 ). Figure H] also shows the mono-mass King 
model that best fits the derived density profile, with the corresponding values of the core 
radius and concentration being r c ~ 9'/7 (with a typical error of ~ ±2") and c = 1.71, 
respectively (hence, the tidal radius is r f ~ 500" ~ 50 r r ). These values are in good agreement 



(11993 . r c = 9"55 and c 



(1996 


, 7> = 9"6 and c = 1.72), 


Trager, Dioreovski & King 


and 


McLaughlin & van der Mare] 


( 


2005 


, r c = 10 / /3 and c = 



1.68), derived from the surface brightness profile, and they confirm that NGC 1904 has not 
yet experienced core collapse. By as suming a distance modulus (m — M) = 15.63 (distance 
d ~ 13.37 kpc, iFerraro et al.lll999bl ). the derived value of r c corresponds to ~ 0.65 pc. By 
summing the luminosities of stars with V < 20 observed within ~ 4", we estimate that the 



extinction-corrected cen tral sur 
in good agreement with lHarris 



McLaug hlin fc van der Marell (12005 



a ce bri ghtness of the cluster is Hy n (Q) — 
Jl996l . hvjq = 16.23), biorgovski I Jl993l . 

fjbv,o = 16.18). Following the procedure described in 



16.20 mag/arcsec 2 , 
Hv,o = 16.15), and 



Djorgovski I (119931 . see also Beccari et al. 2006), we derive logu ~ 
central luminosity density in units of L^/pc 3 (for com parison, logz/ 
Djorgovski Ill993l ; iMcLaughlin fc van der Marell 120051 ). 



3.97, where u n 


is the 


= 4.0 in 


Harris 


1996; 



3. THE BSS POPULATION OF NGC 1904 

3.1. BSS selection 

At UV wavelengths BSS are among the brightest objects in a GC, and RGB stars are 
particularly faint. By combining these advantages with the high-resolution capability of HST, 
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the usual problems associated with photometric blends and crowding in the high density 
central regions of GCs are minimized, and BSS can be most reliably recognized and separated 
from the other populations in the UV CMDs. For these reasons our primary criterion for 
the definition of the BSS sample i s based on the position of stars in the (m 2 is, ^218 — B) 



plane (see also lFerraro et al.l 12004 for a detailed discussion of this issue). In order to avoid 
incompleteness bias and the possible contamination from TO and sub-giant branch stars, 
we have adopted a limiting magnitude m 2 is = 19.5, roughly corresponding to 1 magnitude 
brighter than the cluster TO. The resulting BSS selection box in the UV CMD is shown 
in Figure [5j Once selected in the UV CMD, all the BSS lying in the field in common with 
the optical-HST sample have been used to define the selection box in the (V, B — V) and 
(V, U — V) planes. The limiting magnitude in the V band is V ~ 18.9, and the adopted 
BSS selection boxes in these planes are shown in Figures [3] and [H] (only stars not observed 
in HST-Pointing B are shown in the latter). 

With these criteria we have identified 39 BSS in NGC 1904: 37 in the HST sample 
(32 from HST-Pointing B, and 5 from HST-Pointing A) and 2 in the External sample (r > 
85"), the most distant lying at r ~ 270" (~ 4(5) from the cluster center (see Figure [2]). 
All candidate BSS have been confirmed by visual inspection, evaluating the quality and 
the precision of the PSF fitting. This procedure significantly reduces the possibility of 
introducing spurious objects, such as blends, background galaxies, etc., in the sample. The 
coordinates and magnitudes of all the identified BSS are listed in Table [TJ 

In order to study the radial distribution of BSSs, one needs to compare their number 
counts as a function of radius with those of a population assumed to trace the radial density 
distribution of normal cluster stars. We chose to use HB stars for that purpose, given 
their high luminosities and relatively large number. Thanks to the (essentially blue) HB 
morphology, such a population can be easily selected in all CMDs, and the adopted selection 
boxes, designed to include the bulk of HB and the few post-HB stars, are shown in Figures 
[5HS1 In order to be conservative, a few stars lying within the adopted HB selection boxes 
in the optical bands, but not detected in the UV filters (GALEX-NUV channel and HST- 
F218W filter), have been excluded from the following analysis. However slightly different 
boxes or the inclusion of these stars in the sample have no effects on the results. With these 
criteria we have identified 249 HB stars (197 at r < 85" from the HST sample, and 52 at 
85" < r < r t from the External sample). 
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3.2. BSS mass distribution 



The position of BSS in the CMD can be used to derive a "photometric " estimate of their 
masses through the comparison with theoretical isochrones. We did this in the (V, B — V) 
plane, where 34 BSS (32 from the HST-Pointing B and 2 from the External sample) out of 
the 39 identified in the cluster have been measured. 

A set of iso chrones of appropriate metallicity {Z = 6 x 10~ 4 ) has been extracted from 



the data-base of ICariulo. Pegl'Innocenti fc Castellani I (120031 ) and transformed int o the ob- 
servational plane by adopting a reddening E(B — V) = 0.01 (IFerraro et al.lll999bl ). The 12 
Gyr isochrone nicely reproduces the main cluster population, while the region of the CMD 
populated by the BSS is well spanned by a set of isochrones with ages ranging from 1 to 6 
Gyr (see Figure [7]). Thus, the entire dataset of isochrones available in this age range (stepped 
at 0.5 Gyr) has been used to derive a grid linking the BSS colors and magnitudes to their 
masses. Each BSS has been projected on the closest isochrone and a value of its mass has 
been derived. As shown in the lower panel of Figure BSS masses range from ~ 0.95 to 
~ 1.6M , and both the mean and the median of distribution correspond to 1.2 M & . The 
TO mass turns out to be Mto = O.8M . 



3.3. The BSS radial distribution 

The radial distribution of BSS identified in NGC 1904 has been st udied following 



the same procedur e previously adopted for other clusters (see references in IFerraro 1 12006 



Beccari et al.l 120061 ) . In Figure [8] we compare the BSS cumulative radial distribution to that 
of HB stars. The two distributions are obviously different, with the BSS being more centrally 
concentrated than HB stars. A Kolmogorov-Smirnov test gives a ~ 7 x 10~ 4 probability that 
they are extracted from the same population, i.e. the two populations are different at more 
than 3cr level. 

For a more quantitative analysis, the surveyed area has been divided into 6 concentric 
annuli, the first roughly corresponding to the core radius (r = 10"), and the others chosen 
in order to sample approximately the same fraction of the cluster luminosity out to the tidal 
radius (r t ~ 500"). The luminosity in each annulus has been calculated by integrating the 
surface density profile shown in Figure HI The number of BSS and HB stars (iVess and iVnB, 
respectively), as well as the fraction of sampled luminosity (L samp ) measured in each annulus 
are listed in Tableland have been used to compute the population ratio N BSS /Nn B and the 
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specific frequencies (see lFerraro et al.ll2003l ): 

(^pop 



R 



tot \ 
pop J 



pop 



( L samp/ L samPy 



with pop = BSS, HB. 

The resulting radial trend of Rub over the surveyed area is essentially constant, with a 
value close to unity (see Figure[9]). This is just what expected on the basis of the stellar evolu- 
tion theory, which predicts that the fraction of stars in any post-main sequence evolutionary 



stage is strictly proportional to the fraction of the sampled luminosity (IRenzini fc Fusi Pecci 
19881 ). In contrast the BSS show a completely different radial distribution: as shown in 



Figure [9], the specific frequency Rbss is highly peaked at the cluster center decreases to a 
minimum at r ~ 12 r c and remains approximately constant outwards. The same behavior 
is clearly visible also in Figure [lQl where the population ratio A^bss/^hb is plotted as a 
function of r/r c . 



3.4. Dynamical simulations 



Following the same approach as iMapelli et al. (b004L 2006) and lLanzoni et al. I J2007h 



we ha ve used a Monte-Carlo simulation code (originally developed by lSigurdsson fc Phinney 



19951 ) in order to reproduce the observed radial distribution and to derive some clues about 
the BSS formation mechanisms. Such a code follows the dynamical evolution of N BSS 
within a background cluster, taking into account the effects of both dynamical friction and 
distant encounters. Since stellar collisions are most probable in the central high-density 
regions of the clusters, in the simulations we define COL-BSS those objects with initial 
positions r, ^ r c . Since primordial binaries most likely evolve in isolation if they orbit in 
the cluster outskirts, we identify as MT-BSS those BSS having rj 3> r c . Within these 
defintions, in any given run we assume that a certain fraction of the N simulated BSS is 
made of COL-BSS and the remaining fraction of MT-BSS. The initial positions of the 
two types of BSS are randomly generated within the appropriate radial range (rj r c for 
COL-BSS, and ^> r c for the others) following a flat distribution, according to the fact 
that the number of stars in a King model scales as dN = n(r) dV oc r~ 2 nr 2 dr oc dr. Their 
ini tial velocities are randomly e xtracted from the cluster velocity distribution illustrated 



in 



Sigurdsson fc Phinney I (119951 ). and an additional natal kick is assigned to COL-BSS to 



account for the recoil indu ced by the three-body encounters that trigger the merger an d 
produce the BSS (see, e.g.. ISigurdsson. Davies fc Bolte Ill994l ; iDavies. Benz &: Hills 111994 ). 
Each BSS has characteristic mass M and maximum lifetime ti ast . We follow their dynamical 
evolution in the (fixed) gravitational potential for a time tj (i = 1,N), where each ti is a 
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randomly chosen fraction of ti ast . At the end of the simulation we register the final positions 
of BSS, and we compare their radial distribution with the observed one. The percentage 
of COL- and MT-BSS is changed and the procedure repeated until a reasonable agreement 
between the simulated and the observed distributions is reached. 

For a more detailed discussio n of the procedure an d the ranges of values appropriate for 
the input parameters we refer to iMapelli et al.l (120061 ). Here we only list the assumptions 
made in the present study: 



- the background cluster has been approximated with a multi-mass King model, deter- 
mined as the best fit to t he observed profile^. The cluster central velocity dispersion 
is set to a = 3.9 km s" 1 ( IDubath. Meylan fc Mayor 1119971 ). and, assuming 0.5 M Q as 
the average mass of the cluster stars, the central stellar density is n c = 3 x 10 4 pc -3 
JPrvor fc Mevlan I fl993h : 



- BSS masses have been fixed to M — 1.2 M (see Section I3~2"l) and characteristic lifetimes 
ti as t ranging between 1.5 and 4 Gyr have been considered; 

- COL-BSS have been distributed with initial positions r, < r c and have been given a 
natal kick velocity of 1 x a; 

- initial positions ranging between 5 r c and r t have been considered for MT-BSS in dif- 
ferent runs; 

- in each simulation we have followed the evolution of N = 10, 000 BSS. 

The simulated radial distribution that best reproduces the observed one (with a reduced 
X 2 — 0.1) is shown in Figure [10] and is obtained by assuming that the totality of BSS is 
made of COL-BSS. In the best-fit case the BSS characteristic lifetime is ti ast — 1.5 Gyr, 
but a variation between 1 and 4 Gyr of this parameter still leads to a very good agreement 
(x 2 — 0.2-0.3) with the observations. For the sake of comparison, in Figure [TD] we also show 
the results of the simulations obtained by assuming a percentage of MT-BSS ranging from 
10% to 40% (see lower and upper boundaries of the gray region, respectivelyjE As can be 



4 By adopting the same mass groups as those of IMapelli et all (|2006l ). the resulting value of the King 
dimensionless central potential is Wo = 10 

5 Note that a population of 40% MT-B SS was needed in order to reproduce the bimodal distribution 



Lanzoni et al. 1 120071 ) 



observed in M 3, 47 Tuc and NGC 6752 (|Maperii et al.l 120061) . and 10% was found to be the appropriate 
percentage of MT-BSS in the case of M 5 
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seen, while a population of 10% MT-BSS is still marginally consistent with the observations, 
larger percentages systematically overestimate the BSS population at r £ 5 r c . Increasing 
the BSS mass up to 1.5 M does not change this conclusion. 

By assuming 12 Gyr for the age of NGC 1904, we have used the simulations and the 
dynamical friction timescale (from, e.g., Mapelli et al.l [20061 ) for 1.2 M & stars to estimate 



the radius of avoidance r avoid of the cluster, i.e., the radius within which all these stars are 
expected to have already sunk to the cluster core because of mass segregation processes. We 
find that r avoid ~ 30 r c (i.e., ~ 300"), which corresponds to a significant fraction of the entire 
cluster extension. This evidence is consistent with the fact that the simulated MT-BSS 
appear to be a negligible fraction of the overall BSS population. 



DISCUSSION 



We have studied the brightest portion (m 2 i8 < 19.5) of the BSS population in NGC 1904. 
We have found a total of 39 objects, with a high degree of segregation in the cluster center. 
Approximately 38% of the entire BSS population is found within the cluster core, while only 
~ 13% of HB stars are counted in the same region. This indicates a significant overabundance 
of BSS in the center, as also confirmed by the fact that the BSS specific frequency -Rbss within 
r c is roughly 3 times larger than expected for a normal (non-segregated) population on the 
basis of the sampled light (see Figure [9]). The peak value is in good agreement with what is 



found in the case of M 3, 47 Tuc, NGC 6752 and M 5 (see iFerraro et al.l 12004 ; ISabbi et al. 



2004 ; lLanzoni et al. 1 120071 ). Unlike these clusters, no significant upturn of the distribution 
at large radii has been detected in NGC 1904. 

We emphasize that the absence of an external upturn in the BSS radial distribution is 
not an effect of low statistics. In the case of NGC 6752, where a si milar amount of B SS (34) 
has been detected, the BSS radial distribution is clearly bimodal (ISabbi et al.l 120041 ) . This 
can be seen also in Figure [TT| where the two distributions are directly compared. They nicely 
agree within r ~ 12r c , but the fraction of BSS in NGC 6752 rises again at larger distances 
from the center, despite the smaller number of BSS observed in this cluster compared to 
NGC 1904. 

Extensive dynamical simulations have been used to derive some hints about the BSS 
formation mechanisms. Even if admittedly crude, this approach has been successfully used to 
demonstrate that the external rising branch of the BSS radial distribution observed in M 3, 
47 Tuc, NGC 6752 and M 5 cannot be due to COL-BSS originated in the core and then kicked 
out in the outer regions: hence, a significant fraction (20-40%) of the overall population is 
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required to be made of MT-BSS in these clusters ( IMapelli et al.ll2006l ; lLanzoni et al. 1120070 . 
By using the same simulations to interpret the (flat) BSS radial distribution of NGC 1904, 
we found that only a negligible percentage (0-10%) of MT-BSS is needed. However, we 
emphasize that if a rising peripheral BSS frequency is absent (as in the case of NGC 1904) 
our simple approach cannot distinguish between BSS created by MT (and then segregated 
into the cluster core by the dynamical friction) and COL-BSS created by collisions inside 
the core. 

On the other hand, the negligible fraction of peripheral MT-BSS found in NGC 1904 is 
in agreement with the quite large value of the radius of avoidance estimated for this cluster 
( r avoid — 30 r c ), which indicates that all the heavy stars (binaries) within this radial distance 
have had enough time to sink to the core and are therefore not expected in the cluster 
outskirts. Such a radial distance corresponds to 0.6 r t , i.e., it represents a significant fraction 
of the cluster extension (only 1% of the cluster light is contained between r avoi d and r t ), 
and hence only a small fraction of the massive objects are expected to be unaffected by the 
dynamical friction). In all the other studie d cases, r avoi d is significantly smaller: r avoi d ^ 0.2 r t 
( IMapelli et al.ll2006l ; lLanzoni et al. 1 120071 ). In turn, this suggests that at least a fraction of 
the BSS population that we now observe in the cluster center are primordial binaries which 
have sunk to the core because of the dynamical friction process, and mixed with those that 
formed through stellar collisions. 

Only systematic surveys of physical a nd chemical properties f or a large number of BSS 
in different environments (see examples in lDe Marco et al. 1120051 ; iFerraro et al.ll2006aO can 
definitively identify the formation processes of these stars. 
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Fig. 1. — Map of the combined HST sample. The light solid and dotted lines delimit the FoVs 
of Pointing B and A, respectively. Star positions are plotted with respect to the center of 
gravity C grav derived in Section a(J2000) = 05 h 24 m 11.109, <5(J2000) = -24° 31' 29700. 
The positions of all BSS identified in this sample are marked with heavy dots and the 
concentric annuli used to study their radial distribution (cfr. Table 1) are also shown. The 
inner and outer annuli correspond to r = r c = 10" and r = 85", respectively. 
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Fig. 2. — Map of the External sample. The light solid and dotted lines delimit the WFI and 
the GALEX FoVs, respectively. The two BSS detected in the External sample are marked 
as heavy dots, and the concentric annuli used to study their radial distribution are shown 
as heavy circles. The inner annulus is at 85" and corresponds to the most external one in 
Figure [H The heavy dashed circle marks the tidal radius of the cluster (r t c± 500"). 
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Fig. 3.— (V, B — V) CMDs of the HST {Pointing B) and External samples. The hatched 
regions (V > 20) indicate the stars not used to derive the cluster surface density profile. The 
adopted BSS and HB selection boxes are shown, and all the identified BSS are marked with 
the empty circles. 
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Fig. 4. — Observed surface density profile (dots and error bars) and best-fit King model 
(solid line). The radial profile is in units of number of stars per square arcsec. The dotted 
line indicates the adopted level of the background, and the model characteristic parameters 
(core radius r c , concentration c, dimensionless central potential Wo), as well as the \ 2 value 
of the fit are marked in the figure. The location of the cluster tidal radius is marked by the 
arrow. The lower panel shows the residuals between the observations and the fitted profile 
at each radial coordinate. 



- 21 - 



16 



18 



n 1 r 



20 - 




i 1 1 1 r~ 

HST Sample 
(Pointing B) 




m 2m -B 



Fig. 5. — CMD of the ultraviolet (Pointing B) HST sample. The adopted magnitude limit 
and selection box used for the definition of the BSS population (empty circles) are shown. 
The two solid triangles correspond to BSS-38 and 39 found in the External Sample, with 
UV magnitudes obtained through the GALEX NUV detector. The selection boxes adopted 
for HB and post-HB stars are also shown. 
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Fig. 6. — (V, U — V) CMD of the HST {Pointing A) sample (only stars not observed in 
Pointing B are plotted). The adopted BSS and HB selection boxes are shown, and all the 
identified BSS and HB stars are marked with the empty circles and squares, respectively. 
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Fig. 7. — Upper panel: zoomed (V", B — V) CMD of the BSS region; the 34 BSS measured 
in thi s plane are shown. The set of isochron e s rang ing from 1 to 6 Gyr (stepped by 0.5 Gyr) 
from ICariulo. Degl'Innocenti fc Castellani I (120031 ) data base used to derive BSS masses is 
also shown. Lower panel: derived mass distribution for the BSS shown in the upper panel. 
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log(r/ arcsec) 

Fig. 8. — Cumulative radial distribution of BSS (solid line) and HB (dashed line) stars as a 
function of the projected distance from the cluster center for the combined HST+External 
sample. The location of the cluster tidal radius is marked by the arrow. 
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Fig. 9. — Radial distribution of the BSS (dots) and HB (gray regions) specific frequencies, 
as defined in equation (JT]), and as a function of the radial distance in units of the core radius. 
The vertical size of the gray regions correspond to the error bars. 



-26- 



T 



^ Observations 




10 20 30 

r/r 



Fig. 10. — Radial distribution of the population ratio iV B ss/^HB as a function of r/r c (dots 
with error bars), compared with the simulated distribution (solid line and triangles) obtained 
by assuming 100% of COL-BSS. The results of the simulations obtained by assuming a 
percentage of MT-BSS ranging from 10% to 40% (lower and upper boundaries of the gray 
region, respectively) are also shown. 
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Fig. 11. — Radial distribution of the population ratio A^bss/^hb for NGC 1904 (filled circles) 
and NGC 6752 (open circles) plotted as a function of the radial distance in core radius units. 
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Table 1. The BSS population in NGC1904 



Name 




T> A 
HA 

[degree] 


[degree] 


«^218 


TT 

u 


p. 

-D 


V 


T 
1 


T3QQ 1 
-Dbb-1 


ol 


c\A 07O7i nn 
.U4o l\) ( 1UU 


oa ^okqoi 1 nn 
-z4.ozooyi iuu 


i n 1 1 
iy.ii 


1 Q P.A 
10.04 


1 Q P.P. 

10.00 


~\ Q AK 
10.40 




T3CC O 


O 1 

ol 


.U4f /ozoUU 


-z4.0zb0z / 4UU 


1 o /?c 

lo.oo 


1 O 1 o 

lo.lz 


1 Q Ol 

lo.zl 


i 7 n/i 

1 1 .y4 




T3CC Q 


O 1 

ol 


n /i 7c /i n/i nn 
.U4/04U4UU 


-z4.0zbUU0oUU 


1Q or 

lo.oo 


1 O Q A 

lo.o4 


1 O QO 

lo.oo 


1 O OO 




T3CC A 
-Dbb-4 


O 1 

ol 


n/iooc^/i nnn 
.U4oy04UUU 


oa r;o^i oconn 
-z4.0z01oozUU 


i 7 n /i 

1 1 .y4 


1 7 7C 

1 / . / 


1 7 07 
1 / .0 / 


1 7 70 




T3CC C 

ribb-0 


O 1 

ol 


n /i 71 nnonn 

.U4f lyyzuu 


oa ^o^7n7^nn 

-/4.0Z0 / y I OUU 


1 7 CQ 

1 ( .00 


17/17 

1 / .4 / 


i 7 /I n 
1 f .4U 


1 7 Q C 

1 / .o0 




T3CC a 

ribb-b 


O 1 

ol 


n /i c cooQnn 
.U400zooUU 


o/i c^o^c7^cnn 
-Z4.0Z00 /boUU 


A Q P. A 

lo.04 


1 O O /I 

lo.z4 


1 O OO 

lo.zo 


1 O 1 o 

lo. iZ 




T3CC 7 

r3bb- ( 


O 1 

ol 


n /i /i on^ 1 nn 
.U44zyblUU 


-z4.0zbobZZUU 


1 n 1 q 

iy.io 


1 O 77 

lo. / / 


1 O 7£ 

lo. / 


1 o c c 
lo.OO 




T3CC C 

rJbb-0 


Q 1 

ol 


n/i q tcnfi7nn 
.U4oOUb / UU 


O/i co/ioKfjonn 
-z4.0z4zbboUU 


i n i q 

iy.io 


1 Q C C 

lo.OO 


IOC/I 

lo.04 


i q Qn 
lo.oy 




T3CC n 

rJbb-y 


O 1 

ol 


n /i 1 c c^ 1 nn 
.U410001UU 


oa c;ofin70/inn 
-z4.ozby l Z4UU 


1 n o c 
ly.oO 


lo.OO 


1 O 0£ 

lo.ob 


1 q /i n 

io.4y 




T3QQ 1 n 


ol 


n/i coo7Qnn 
.U400Z I oUU 


o/i ^oc;n£:o7nn 
-Z4.0Z0U0Z ( UU 


i s nf; 
lo.Uo 


17 QC 
1 1 .OO 


17 17 

It A i 


1 of; 
Ib.yb 




T2CO 1 1 

-Bbb-11 


Q1 

ol 


n/i K/ifi77nn 
.U4040 1 t\j\j 


O/i co ci /i 7nnn 

-z4.ozoi4 1 yuu 


1 7 C A 

1 / .04 


1 7 /I A 

1 / .44 


1 7 C£ 

1 1 .0b 


17/11 

1 / .41 




rJbb-lz 


01 

ol 


.U4<UoooUU 


O/i co/iovcnnn 


i o on 
lo.oU 


1 o on 
lo.zU 


i o m 
lo.Ul 


i 7 on 
1 / .oU 




-Dbb-lO 


Ol 

ol 


c\ap^aq a nn 
.U4b04o4UU 


O/i co/i/ioovnn 
-z4.0Z44oo / UU 


i n /i q 

iy.4o 


1 O 77 

lo. I I 


1 O ^7 

lo.b / 


1 O QO 

lo.oo 




T3CC 1 A 

-bob- 14 


01 

ol 


n a c i o 1 onn 
.U401ZloUU 


oa co/17/ioonn 
-Z4.0Z4 / 4ooUU 


i n oo 

iy.zo 


1 O QQ 

lo.oo 


lo.OO 


1 O Ol 

lo.zl 




T3QQ 1 c 
ribb-10 


Q1 

ol 


n/i Kooqqnn 
.U40ooooUU 


OA KO AO10.RC\C\ 
-Z4.0Z4Z I OOUU 


1 fi 07 

lo.y i 


~\Q AO 
10. 4Z 


1 q Qn 
lo.oU 


1 Q 1 Q 

lo. lo 




rJbb-lb 


01 

ol 


r*AP~\ p. A enn 
.U4blb40UU 


O/i cooc^7Qnn 
-Z4.0ZZ00 / oUU 


1 O A 1 

10.41 


1 7 70 
Li.il 


I 7 £Q 

I I .bo 


1 7 /I O 

1 r .4z 




T3CC 1 7 


C 1 

ol 


n/i a qi Q7nn 
.U44olo / UU 


OA KOQQfl/IOnn 

-z4.ozooU4zUU 


lo.oo 


1 Q /I P. 

lo.4b 


10/1/1 

10.44 


1 O 07 

lo.z 1 




-Dbb-lo 


C 1 

ol 


.U4< 10 /OUU 


oa coi noo7nn 
-z4.0Ziyoz / UU 


i o /I n 
lo.4y 


1 O /I c 

1O.40 


10/11 

10.41 


1 O QO 

lo.oz 




T3CC in 

-Dbb-ly 


ol 


n /i o /i i oi nn 
.U4o41olUU 


oa KOQO/iT'onn 
-z4.ozoz4 1 zUU 


1 7 OO 
1 / .00 


lo.lb 


I 7 £0 

I I .bo 


1 7 C C 

1 / .00 




Dbb-Z(J 


01 

ol 


n /i £££ rnnn 


oa com nonnn 
-z4.0ZUiyzUUU 


i n /i n 

iy.4y 


1 O 71 

lo. 1 1 


1 o no 

lo.yz 


1 o ^n 
lo.OU 




ncc 01 
Dbb-Zl 


01 

ol 


c\A a nm /i nn 

.U44yy i4uu 


o/i com ovenn 
-z4.0ZUlZ 1 OUU 


i o no 
lo.yo 


1 O /I c 

1O.40 


1 O A A 

10.44 


1 O Ol 

lo.zl 




T2CO OO 


ol 


.U4010 /oUU 


O/i ci no /i ci nn 
-z4.0iyz4olUU 


i a aci 

iy.4y 


1 o no 

lo.yz 


1 n c\p 

iy.ub 


1 O 7/1 

lo. /4 




.Dbb-ZO 


Ol 

ol 


n a noo^ 1 nn 
.U4yoZOlUU 


oa coi £01 nnn 

-z4.ozioziyuu 


i o no 
lo.Uo 




i o n 1 
lo.Ul 


1 7 nn 

i ( .yy 




ncc oa 
rJbb-z4 


ol 


n /i ni c ennn 

.u4yiooyuu 


oa con^onenn 
-z4.0ZUozyt)UU 


1 O /I 1 

10.41 




i 7 no 

i ( .yo 


1 7 07 
1 1 .0 / 




DCQ oir 

rJbb-zO 


C 1 

ol 


n/i 70/1 /i Knn 
.U4/Z440UU 


oa ci 7ncot?nn 
-z4.ol /UozoUU 


i n a i 

iy.4i 


1 Q (17 

lo.y / 


1 n 1 1 

iy.ii 


1 Q QO 

lo.oz 




pec op 
ribb-zb 


1 

ol 


nri cnoonn 

.uoioyzzuu 


O/i coqi /ififinn 
-z4.ozol4ooUU 


1 Q R7 

lo.b / 




1 O QQ 

lo.oo 


1 O 01 

lo.zl 




T2CC Q7 
.Dbb-Z 1 


CI 

ol 


n cn a *?p 1 nn 
.U0U4/01UU 


O/i ci 71 nvnnn 
-Z4.01 1 1U / UUU 


1 n n 1 

iy.ui 


lo (?r 

lo.OO 


1 Q C A 

1O.04 


i o Qn 
lo.oy 




T3CC OQ 
-Dbb-ZO 


ol 


.UbU fb rfUUU 


O/i conncQnnn 
-z4.ozUyooyUU 


1 o no 
lo.yz 




i o cn 

lo.oy 


1 Q /I C 

1O.40 




BSS-29 


81 


.068117100 


-24.517558600 


19.39 




18.91 


18.60 




BSS-30 


81 


.043233100 


-24.533852000 


17.34 




16.75 


16.65 




BSS-31 


81 


.044917900 


-24.540315000 


18.34 




17.48 


17.22 




BSS-32 


81 


.038520800 


-24.540891600 


17.77 


17.46 


17.30 


17.20 




BSS-33 


81 


.045196533 


-24.515874863 




17.91 




17.84 




BSS-34 


81 


.032196045 


-24.512256622 




18.75 




18.41 




BSS-35 


81 


.037574768 


-24.525493622 




18.87 




18.58 




BSS-36 


81 


.041069031 


-24.518457413 




18.80 




18.64 




BSS-37 


81 


.045227051 


-24.517648697 




18.86 




18.68 




BSS-38 


81 


.056510925 


-24.449676514 


19.04t 


18.78 


18.56 


18.35 


18.08 



-29- 



Table 1 — Continued 
BSS-39 81.058883667 -24.555763245 19.39 f 19.14 19.05 18.73 18.34 



Note. - - t Note that, while the header of the column referes to HST-F218W 
magnitudes, those of BSS-38 and -39 have been obtained with the NUV channel 
of GALEX and transformed to the m 2 is scale as described in Section 2.2. 



n 


r e 


^Vbss 


iV HB 


r samp / r samp 
^ 1 ^tot 





10 


15 


34 


0.14 


10 


20 


10 


45 


0.18 


20 


40 


7 


62 


0.22 


40 


85 


5 


56 


0.23 


85 


150 


1 


34 


0.13 


150 


500 


1 


18 


0.10 



Table 2: Number of BSS and HB stars, and fraction of luminosity sampled in the 6 concentric 
annuli used to study the BSS radial distribution of NGC 1904 (r^ and r e correspond to the 
internal and external radius of each considered annulus, in arcsec). 



